LOCAL GAUSSIAN FLUCTUATIONS IN THE AIRY AND 
DISCRETE PNG PROCESSES 



JONAS HAGG 

Abstract. We prove that the Airy process, A(t), locally fluctuates like 
a Brownian motion. 

In the same spirit we also show that, in a certain scaling limit, the so 
called Discrete polynuclear growth process (PNG) behaves like a Brow- 
nian motion. 



1. Introduction 

1.1. The Airy process. The central object of study in this paper is the 
local behavior of the Airy process, t — > A(t), t G R, p3]. The Airy process 
is a one dimensional process with continuous paths, [6J, [13]. The interest 
in this process is mainly due to the fact that it is the limit of a number of 
processes appearing in the random matrix literature. One example is the 
top curve in Dyson's Brownian motion, see [3], which, when appropriately 
rescaled, converges to the Airy process, see for instance [2] and [7J. Another 
example is the boundary of the north polar region in the Aztec diamond, 
see [1], [5] and [8], a discrete process also converging to the Airy process, 
[Sj. A third example, the Discrete polynuclear growth model (PNG), [7], [9], 
will be described in some detail in section 11.31 where we also state a theorem 
about its local (in a certain sense) fluctuations. 

A precise definition of A(t) goes as follows: 

The extended Airy kernel, [2], [10], [13], is defined by 



(1.1) A s>t (x,y) 



J °° e- z ^Ai(x + z)Ai(y + z)dz if s > t 
- / M e^-^Aiix + z)Ai(y + z) dz if s < t, 



where Ai is the Airy function. A StS (x,y) is easily seen to be the ordinary 
Airy kernel, [15]. Given £i, . . . , £ m £ R and ti < . . . < t m in R we define / 
on {tx, . . . , t m } x R by 

f(U,x) = X(^,oc)(x). 

It is shown in [7J that 

f 1 /\s,x)A s>t (x,y)f 1 / 2 (t,y) 

is the integral kernel of a trace class operator on L 2 ({ti, . . . , t m } x R) where 
we have counting measure on {t\, . . . , t m } and Lebesgue measure on R. The 
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Airy process, t — > A(t), is the stationary stochastic process with finite di- 
mensional distributions given by 

F[A(t 1 )<Z 1 ,...,A(t m )<U]=det(l-f 1 / 2 Af 1 / 2 ) 

The determinant in the right hand side is a so called Fredholm determinant. 

Our main theorem states that if we condition the Airy process to be at 
some given point at time t\ it will then behave, on a local scale, like a 
Brownian motion. 

Theorem 1.1. Let e > be small, t\ G M. and tj = t{-\ + s«e, 2 < i < m, 
where s 2 , • • • , s m > 0. Also, let p\ G M. and define the sets At, i = 2, . . . , m, 
by 

Ai = [x G R|pi + aj^/e < £ < Pi + biVe} 
where ai,bi are given real numbers. It holds that 



P [A(t 2 ) G A 2 , . . . , A(t m ) G A m \ A(h) = Pl ] 



&2 rbm 

dx 2 ■ J dx m 

a,2 J am 



2 m 
E 2 



\/4vrs2 



n 

i=3 



e 4s i 



+ E 



TTSi 



where 



l^l^v^loge 1 JJOi - a,i)C pitS2 ^ s 



i=2 



Figure 1 describes the setup in the theorem. 
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Figure 1. Conditioned that A(t\) = p\ Theorem ll.il gives 
the approximate probability for the process to move through 
the sets A{. Note that — t j ~ e and \ A^\ ~ ^fe. 



Remark 1. A couple of previous results about the Airy process are the 
following: 

In [13] it is shown that 

Yar(A(t) - .4(0)) = 2t + 0{t 2 ) 

as t — » 0. 
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In p] , see also [16] , the long distance covariance asymptotics for the Airy 
process is calculated to be 

E [A(t)A(0)] - E [Ait)] E [A(0)] = r 2 + 0{t~ 4 ) 

as t — > oo. This proves that A(t) is not a Markov process since this would 
imply exponential decay. 

Remark 2. Given Theorem 11.11 it is natural to ask the corresponding 
question about processes converging to the Airy process. Theorem 11.31 in 
section [L3l below provides such a result for the Discrete polynuclear growth 
process. 

1.2. The extended Airy point process. We now present another con- 
struction, [TJ, of the Airy process that will help us analyzing its local behav- 
ior. 

Let m G Z + be arbitrary and t\ <t%... < t m be points in R which we 
shall think of as times. Define 

£ = l !l UM( 2 U.-Ul im . 

We shall refer to Rj. as time line tj. We define X to be the space of all 
locally finite countable configurations of points (or particles) in E. Locally 
finite means that, if x = (x\,X2, ■ ■ ■) G X then, for any bounded set C C E, 
it holds that #(C D x) < oo. Here #5 represents the number of points in 
the set B. One can construct a a- algebra on X from the cylinder sets: Let 
B C E be any bounded Borel set and n > 0. Define 

to be a cylinder set and X to be the minimal er-algebra that contains all 
cylinder sets. One can now define probability measures on the space (X, £). 
The so called extended Airy point process is an example of such a measure 
and it will be described below. 

For the sake of convenience, we will often denote the extended Airy kernel 
by A(x,y) instead of At it f.(x,y) when it is clear that x G Mf i and y G 
Let Z\, ,Zk be points in E. The fc-point correlation function is defined by 

(1.2) R(z u ...,z k ) = det [A(zi, «i)]f J=1 • 

It is possible to show that these correlation functions determine a probability 
measure on (A, S), the extended Airy point process, by demanding that the 
following identity holds, [LI] : 



(1.3) E 



i^Bi-hy. 



/ , R(zi, . . . , z k ) dz. 



Here B 1: ...,B n are disjoint Borel subsets of E and A?j G 1 < i < are 
such that k\ + . . . + k n = k. 

It is possible to show that, at each time line R^, there is almost surely a 
largest particle, A(ij), and 

(1.4) (A(tx),...,A(t m )) = (.A(ii),...,.4(t m )) 
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in distribution, [7]. It is through this representation that we are able to 
show that the Airy process behaves locally as a Brownian motion. 

1.3. Discrete polynuclear growth. The second object of interest in this 
paper is the the so called Discrete polynuclear growth model (PNG), [7], [9]. 
It is defined by 

(1.5) h(x, t + 1) = max (h(x - 1, t), h(x, t), h(x + 1, t)) + u(x, t + 1) 

where x G Z, t G N, h(x,0) = Vx G Z and w(x,i + 1) = if \x\ > t 
or if t — x is even, otherwise u(x, t + 1) are independent geometric random 
variables with 

(1.6) F[u)(x,t + 1) = m] = (1 - q)q m < q < 1. 

It is convenient to extend the process to all x G K by setting /i(x, t) = 
/i([a;J,t). A description of this process using words and pictures goes as 
follows: 

At time t = 1 a block of width one and height u(0, 1) appears over the 
interval [0, 1). This block then grows sideways one unit in both directions 
and at time t = 2 two blocks of width one and heights a; (—1,2), a; (1,2) 
are placed on top of it over the intervals [—1,0) and [1,2) respectively. 
These blocks now grow one unit in each direction disregarding overlaps. 
At time t = 3 three new blocks are placed over [—2,-1), [0,1) and [2,3). 
This procedure goes on producing at each time the curve h(x, t) that can be 
thought of as a growing interface. Figure 2 shows a realization for t = 1, 2, 3. 



-1 



1 

t = 1 



1 

t = 2 



-2 -1 



1 

t = 3 



Figure 2. A sample of the discrete PNG process for t = 
1,2,3. The shaded blocks represent the growth due to the 
random variables uj(x,t). 

The process h is closely connected to a growth model, G(M,N), stud- 
ied in [6]. Let w(i,j), G Z^_, be independent random variables with 



GAUSSIAN FLUCTUATIONS IN THE AIRY PROCESS 5 

distribution given by (II. 6|) . Define 

G(M, N) = max V w(i,j) 

(ij')G7T 

where the maximum is taken over all up/right paths from (1, 1) to (M,N). 
One can think of G(M, N) as a point to point last-passage time and 

G„i(N) = max G(N + K,N- K) 

\K\<N 

as a point to line last-passage time. In [7J it is shown that 

G(i,j) =h{i-j,i + j-l). 

The definition of G p i therefore inspires the study of K — > h(2K,2N — 1), 
that is, the height curve at even sites at time 2N — 1. 

In [7J the rescaled process, t — > H]y(t), t G R, is, for appropriate t, defined 

by 

dN l '*H N (t) = h ( 2 l -±^d- l N 2 /% 2N-l)--^ 



and for the rest of M by the use of linear interpolation. The constant d is 
given by 

d _ U/g^l+y/g) 1 / 3 

i- V? 

The main result about Hjq{t) in [7] is the following theorem: 

Theorem 1.2 (Johansson). Let A(t) be the Airy process defined by its finite 
dimensional distributions and T be an arbitrary positive number. There is 
a continuous version of A(t) and 

H N (t) -» A(t) - t 2 

as N — ► oo m £/ie weak* -topology of probability measures on C(—T,T). 

In particular this theorem shows that the fluctuations of h are of order 
iV 1 / 3 and that non-trivial correlations in the transversal direction show up 
when looking at times ti where ti+\ — ti ~ jV 2//3 . 

Motivated by Theorems 1 1 . 1 1 and 1 1 . 2 1 one could guess that h, on a time scale 
of order iV 7 , < 7 < 2/3, behaves like a Brownian motion. The theorem 
below shows that this is indeed the case. 

Given some m £ Z + set 

= 1 + ^^-1^2/3 

*Q +1 = ^ + i+v^^i ^ i = 1 m _ 1 

where < 7 < 1 and n, Sj > are real numbers such that K{ E Z. Define 



= JVL^ + ^^1/3 G 
1 - Vq 
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where ip is any real number such that J\ G 
Theorem 1.3. Define the sets A{, i = 2, . . . , m, by 

At = [j e Z + \j = Jx + XidN^ 2 , on < Xi < 6i} 

where 04, bi are given real numbers. There exists c > such that 

F[h(2K 2 ,2N -1) G A2, h(2K m , 2N — 1) eA m 

\h(2K 1 ,2N- 1) = Ji] 



62 



1 2__ T-r 1 



dx 2 --- dx m e 4s 2TT— =e 4s ; + £ 



where 



\E\ ^N-'H^-a^C^ 



,S2,.. 

i=2 



2. Proof of theorem 11.11 

The connection (jl.4p shows that we can prove the theorem by studying 
the largest particle in the extended Airy point process at times t\,...,t m . 

The appearance of C in formulae below should be interpreted as follows: 
There exists a positive constant which may depend on pi, Sj, i = 2, . . . , m, 
validating the inequality to the left when inserted instead of C. Other error 
terms will typically also depend on pi, Sj. 

Set Ji = [pi - 61, pi] C R tl and Ji = \pi - ^5i,pi] C R fi , 2 < i < m, 
where 5i > and pi = p%-\ + y%\/e, y% € M. We also set h = (pi, 00), 
i = 1, . . . , m. 

We will show that 

f y 1 P[#Jl > l,...,#J m > l,#Jl = ...= J m =0] 

1 ' ^ 5a,... .S-.0+ «fe • • • *m P > 1, #h = 0] 

1 _ J£2__ »m 

e 4s 2 "' + O(Veloge), 



implying Theorem ll.il 

The first step is to show that the probabilities in the numerator and 
denominator above can be approximated by appropriate expected values. 

For k,n £ Z + we shall use the common notation 



n [k] = n(n - 1) • • • (n - k+ 1). 
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Let J be an interval on some time line and xa be the indicator function for 
the event A. Since 

( k-1 -#J=k>2 
# J -X{#j>i} - | . #J = o,l 

#JM=#J(#J-1) = If-V \t J J = = k ^ 2 

it holds that 

(2.2) 0<#J-x { #j>i}<#J [2] . 
This together with the following facts will be useful: 

(2.3) P[#J 1 >l,...,#J m >l,#J 1 = 0] 

- P [# J x > 1, . . . , # J m > 1, #/ x = . . . = #I m = 0] = 

= P[#Ji > ,#J m > = 0,(#J a = ... = #I m = 0) c ] 

= P [#Jl > 1, . . . , # J m > 1, #/l = 0, lj£ 2 {#/ 4 / 0}] 
< ^ P [# J X > 1, . . . , # J m > 1, 

We now express the probabilities in terms of expected values. If we set 

(2-4) T(J t ) = #J l - X{#Ji > 1} . 

then 

P[#Ji > l,...,#J m > = 0] 

= E [(#Jx - T(Jx)) • • • (#J m - T{J m )) ■ X{#i 1= o}} 
= E [(#Jx ...#J m + U(J u ..., Jm)) • X{#/ 1= 0}] 

where C/ is defined by the last equality. In view of (|2.2p and (|1.3p we get, 
for example, 

E pVi) • #J 2 • • • #J m ] < E [# j| 2] • # J 2 • • • #J m 

= / R(x 1 ,x 2 , . . . ,x m+ i)dx = 0(5j ■ 5 2 - ■ ■ 5 m ). 

Since U{J\, . . . , J m ) is a sum of terms like this one (at least one T(Ji)) 
we see that 

x 1 ? n n+ A 1 A E [^( J l' • • • ' J ™) " ^{#^1=0}] = °- 
fll,...,fl m ^U T oi • • • o m 
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Repetition of this argument shows together with (|2.3H that 
1 



lim 



5i,...,<5m->0+£i ■■■S„ 



1 [# Ji > 1, . . . , # J m >l,#h = ... = I m = 0] 



lim 



1 



/ m 

+ o ( lim 



-E [#Ji • • • #J m ■ x { / 1= o}] 



v i=2 



5 1 ,..., 1 5 m -^0+ o"i • • • S r 



-E [# Ji • • • # J m • X{#i^#h] 



and also that 



lim lp [# Ji > 1, #/! = 0] = lim 1e [# Ji • X{ h =o}} ■ 
5^0+ Oi 5^0+ 0i 



Later it will be shown that 

1 



(2.5) 



lim 



8 1 ,...,8 m ^0+ Si---6n 



■E [#Ji • • • #J m • = O(^loge) 



but let us first be constructive. 
We want to show that 



(2.6) 



lim 



1 



s u ...,s m -*o+ 8i---5 T 



-E [# Ji • • • # J m ■ X{Jl= o}] 



lim -E [# Ji • X {h=oy] 
5i^0+ Oi L 



»I _ Vm 



v / (47r) m - 1 s 2 ---s r 



To start with we need to find a representation of the left hand side of (|2.6 
that is suitable for analysis. 



E[#Ji---#J m -X{/ 1=0 }] =E 

= E 

= E 



# Ji • • • # J m • lim e~ x * h 



A— >oo 



(e- A ~ l) k rW 



# Ji • • • # J m • hm V v " L , ~ y I[ 



fc=0 



E 

£;=0 



(-i) fc 



fc=0 



fc! 



-E 



# Jl • • • # J m • #/[ 



In the second equality we have used the formula 



(2.7) 



„An 



E 

fc=0 



(e A - 1)* 



-n 
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In the fourth equality we take the sum out of the expectation. By Fubini's 
theorem we are allowed to do this since 



E 



#ji • • • #j m • ^2 



k=0 



< E 



00 #I k 



k=0 



E 



# h ■ ■ ■ #J m ■ e# h ] < E • • • #Ji] 1/2 E [e 2 #^ 

In fact E [z^ 1 ] is an entire function in z, |14j . 

Another technical issue we need to deal with is to prove that 



1/2 



< 00. 



lim 



(-1)' 



5i,...,6 m -*0+ Sl - ■ S m ^ kl 



E 



E 

k=0 

00 

E 

k=0 



lim 



#Ji • • • #J m • #I t 
1 



fc! <f 1 ,...,5 m -+o+ St ■ ■ ■ S v 



-E 



#Jl • • • # J m ■ #I{ 



[k] 



kl J ik 



(\/i) m R(pi,...,p m ,xi,...,x k )dx. 



Please recall definition (|1.2p and note that the second equality is immediate 
from (TO) . Define G k (zi, z m ), Zi G R ti , by 



(2.8) G k ( Zl ,...,z m ) = 
The identity sought for is 



-1) 



kl 



R(z 1 ,...,z m ,xi,...,x k )dx. 



(2.9) lim 

8 u ...,6 m ->0 + d l ■ ■ ■ ° 



00 

^7-E 



fc=0 



G k (zi,. . . ,z m )dz 



JlX---X J„ 



= E(^) m - 1 G' fc (pi,...,p m ). 

fc=0 

This will hold if for some neighbourhood f2 of (p%, . . . ,p m ) there exist con- 
stants C k > such that 

\G k {zi, . . .,z m )\ < C k 

if (z\, . . . , Zfy%) € f2 and 

00 

y^Cfc < 00. 

fc=0 

That this is indeed the case follows from calculations similar to the ones 
appearing in the proof of Lemma [2 . 2 1 which is given at the end of this section. 
The following lemma can be found in |llj : 



Lemma 2.1. Let a > 0, then 



e az Ai(x + z)M(y + z) dz 



v 7 ! 



-_e 4a 



ira 
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In this section we call this function cj) a (x,y) or simply 4>(x,y) when it is 
clear what a is. From Lemma 12.11 and the definition of the Airy kernel it 
follows that, for s < t 

A Stt (x,y)= / e z ^~^Ai(x + z)Ai(y + z)dz-^ t _ s (x,y) 
Jo 

-■ A s ,t(x,y) - 4>t- s {x,y). 
For s > t it is convenient to set A s j(x,y) = A s j(x,y). 

Lemma 2.2. Suppose that 1 < v < m, v £ Z. Then, for some C depending 
on pi,... ,p m , 

(\/e) m_1 / R(pi,---,Pm,x 1 ,...,x k )dx 

(2.10) = (V~e) m ~ 1 <f>(pi,P2)<P(p2,P3) ■ ■ ■ <f>(Pm-l,Pm) 

x J^R( Pl ,x 1 ,...,x k )dx + ^0 ((Ck)^ 

Furthermore, ifv>2 then 

(Ve)" 1 - 1 [ dx / dyR( Pl ,..., Pm ,x,y) 
Jh Ji v 

(2.11) = (v^)™" 1 0(P1,P2)^(P2,P3) • • • 4>(Pm-l,P m ) 



(^J^R(p 1 ,x l ,x 2 )dx + R(p u x)6xj +0(v^loge). 



From (pnUj) we now get ([2T5]) . 
We turn now to (|2.5|) . Clearly 



E [# Ji ■ ■ ■ # J m • < E [# Ji • • • # J™ • (#/< - #/i 

= E 

We now obtain (12.51) since 



#Jl • • • #Jm ■ (#if + + #h + #/l " 2#/i#/i 



(\/e) m 1 ^ ^ #(Pl , • • • , Pm, x, y) dxdy + J R{p x ,...,p m ,x,y) dxdy 
(2.12) + / R(pi,...,p m ,x)dx+ / R(p 1 ,...,p m ,x)dx 

-2 / R(pi, ■ ■ ■ ,Pmi x, y) dxdy ) =0(Veloge) 

by Lemma 12.21 

To get (|2.ip we need one more result, namely that 



(2.13) lim 1e[#J iX{# / 1= o}] > 0. 

01— >0+ 0\ 
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Let F 2 (s) be the Tracy- Widom distribution function corresponding to the 
largest eigenvalue in the GUE, [15]. Then 



1 . ^ (—1) f f 
(2.14) = lim ^V^-p- / dxQ \ d k x det(A(xi, Xj)) 0<i:j<k 

5 ^ 0+ <*i t^o hl J * Jl i 

= Yl Hp / k ^t(A(x i ,x j )) < iij < k d k x = i^(pi) 
fe=0 ' 

where in the last row xo = p\- The last equality can be obtained by dif- 
ferentiating the corresponding equality for the distribution function F 2 (t), 
|15j . we omit the details here. The first equality has been shown above and 
the second is a special case of (|2.9p . Since i^(s) > for all s£l, see [15] . 
we obtain (|2. 13|) . 

What is still left is to prove Lemma 12.21 
Proof of Lemma \2.2\ We start with (|2.10p . For < r < m — 1 and k > 1 
define D r (k) by 



D r (k) = {V^Y4 l {Pl,P2)<P{P2,P3) ■ ■ ■ <P{Pr,Pr+l) / dx 

J Ik 

A(Pr+i,Pi) VeA(p r+1 ,p r+2 ) ... y/eA{p r+1 ,p m ) A(p r+1 ,Xj) 

A{p m ,pi) ^A(p m ,p r+2 ) ... y/lA[p m ; Pm 
A(x i} pi) ^A(xi,p r+2 ) ... ^feA{xi,p m ) A(xi,Xj) 

In the determinant 1 < i,j < k and for r = we set the empty product 
in front of the integral to 1. Please note that Do(k) is equal to the left 
hand side in (|2.10p . We let D r (k) be almost the same as D r (k). The only 
difference is that we put in A(p r+ i,p r+2 ) in position (1,2) in the matrix 
instead of A(p r+ i,p r+2 ). By using induction we shall now prove that 

(2.15) Do(fc) = D r (k) + ^~eO {{Ck)^ 



for < r < m — 1. Clearly (|2.15p holds if r = 0. Suppose now that (|2.15p 
holds for some r such that < r < m — 2. By expanding the determinant 
in D r (k) along the first row we see that 

(2.16) D r (k) = D r+1 (k) + D r (k). 

What has to be proved is hence that 

D r (k) = yfeO i{Ck)^ 



To do this, Hadamard's inequality will come in handy but before we re- 
call this inequality we present a lemma which will be frequently used from 
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now on. The proof is readily obtained from Lemma 12.11 and the standard 
estimates, see [12] . 

|Ai(x)| < C Af e- 2 N 3/2 / 3 
|Ai'(*)| < C M ^/W\e- 2 ^ i/2 ^ 

that hold for x > —M. 

Lemma 2.3. Suppose that s < t and M > 0. For x,y > —M and any 
A > it holds that 

\A t , s (x,y)\<C M ,xe- Xix+y) 

A t ,.(x, y) = A ttt {x, y) + 0(t - s) e~ x{x+ ^ 

A s>t (x, y) = A t , t (x, y) - (1 + <D(t - a)) 

+ 0(t- s)e~^ x+y) . 

The errors depend only on M and A. Moreover, 

\A SjS (x + a,y)- A s , s (x,y)\ <aC M , A e- A ^ 

for all a > 0. 

Let B = (bi t j)i<i t j< n , bi,j G M be a matrix. Hadamard's inequality states 
that 

1/2 



(2.17) |detS| < IIE 6 



n n 

2 

'ji 

i=l j=l 



Below we find upper bounds for the equivalent to Yl^=i ^fi ^ n the ma trix 
appearing in D r (k). 
Column 1: 

m k 

A 2 (p r+1 , Pl ) + Y,A 2 (x j , Pl ) <C(k + m) 

j=r+l j=l 

Column 2: 

(m k 
A 2 (p r+1 ,p r+2 ) + ^ A 2 (pj,p r+2 ) + ^2A 2 (xj,p r+2 ) 
j=r+2 j=l 

C(k + m) if v > r + 2 

Cm + CJ2 k j =i{A{x j ,p r +2) ~ 4>{xj,p r+ 2)) 2 if v < r + 2 

Columns 3, . . . , m — r (r + 3 < i < m): 

(m k \ 

j=r+l j=l J 



m 
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A 2 ( Pj ,x i ) + Y J A2 ( x j^i) 

j=r+l 3=1 



< J Ep+i [Afaxi) - 0(p is Xi) J + Cfce" 2 *' if u > r + 2 
( C(k + m)e~ 2xi if v < r + 2 

Next we multiply everything together, take the square root and then inte- 
grate. Assume that v < r + 2. 



C(k + m)e [ C + Cj^fo.p,.^) - ^(x j} p r+2 )) s 



x (C(m + k)) m - r - 2 (C(k + m)) fc e- 2 ( Xl+ - +x *> 



< v^(Cfc)^ ^ e -(^+-+^) f i + ^(1 + 



The case t> > r + 2 can be treated similarly. 
To obtain (|2.10j) it remains to show that 



1/2 



dx 



(xj,p r+2 )) dx 



< ^(ckY 



det 



Pi) ^(p 
A{xi,pi) A(xi,Xj) 



dx 



det 



l<i,j<& 
A(Xj,Pl) A(Xi,Xj 



dx + ^0 ((Ck) 



2 



i<t,i<* 

This is quite easily achieved using Hadamard's inequality and Lemma | 
We do not present the details here but instead go on to prove (|2.1ip . 

The first part of the proof will be similar to the proof of (|2.10p and the 
second part is an application of Lemma 12.41 below. 

Let D r {2) and D r (2) be as defined above with the exception that the 
variables x\ and x 2 are now integrated over I\ and I v respectively. By 
construction Dq(2) equals the left hand side in (|2.1ip . If we can show that 

(2.18) D r (2) = 0( v / i) 

then by the same argument as above 

D (2) = L> m _ 1 (2)+0 vv / i). 
To see this we shall only need the trivial fact that 

n n 

\detB\ < Yl^2\bji\ 

i=l j=l 
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where as before B is a real n x n matrix. Define B as the (m + 2 — r) x 
(m + 2 — r) matrix appearing in D r (2). We now estimate the column sums 

n 

Bi := ^ \bji\. 

3=1 

Column 1: 

B 1 = \A( Pr+1 , Pl )\ + ... + \A( Pm , Pl )\ + \A{x uPl )\ + \A{x 2 , P i)\ < Cm 
Column 2: 

B 2 = Ve \ \A(p r+1 , Pr+2 )\ + \A(p r+2 , Pr+ 2)\ + ■■■ + \A(p m ,p r+2 )\ 

+ \A(x 1 ,p r+2 )\ + \A(x 2 , Pr+2 ) 

< (Cm+ \A(xi, Pt+2 )\ + \A(x 2 , Pr+2 )\) 
Middle columns (if any) (r + 3 < i < m): 
B i = V~e (\A( Pr+ i, Pi )\ + ... + \A( Pm , Pi )\ + \A( Xl , Pi )\ + \A{x 2 , Pi )\) < Cm 
Last two columns: 

B m -r+l = \ A( Pr+1 ,Xi)\ + . . . + \ A( Pm ,Xi)\ 

+ \A(xi,X!)\ + \A(x 2 ,xi)\ < Cme~ xi 
B m _ r+2 = \A(p r+1 ,x 2 )\ + . . . + \A(p m ,x 2 )\ 

v-l 

+ \A( Xl ,x 2 )\ + \A(x 2 ,x 2 )\ < Ce~ X2 +</>(x 1 ,x 2 ) + Yl ^(Pfe^a) 

fc=r+l 

Consider the estimates above for £> 2 and B m _ r+2 . The function A(x2, Pr+2) 
will contain a ^-function if and only if v < r + 2, but in this case the sum 

u-l 

k=r+l 

is empty. This means that we do not get terms like 

4>(x2,Pr+2)4>(Pk,X2) 

in the product B 2 B m - r+2 . Given this observation it is easy to see that 
/ B 2 B m ^ r+ iB m ^ r+ 2 dx = 0{\fe) 

J IlXlv 

and this proves (12.181) . 
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The second part of the proof consists of showing that 
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(2.19) 



det 



hy.lv 



A(p 

Pi) A(p 

in i m-1%2) 

A(xx,pi) A(x 1 ,x 1 ) A(xi,x 2 ) 

A(x2,pi) A(x 2 ,xi) A(x 2 ,X 2 ) 



dx 



= / R(p\,xi,x 2 ) dx+ / R(pi,x) dx + 0(y/eloge). 

Jq Jh 

The left hand side is equal to 



det 



hxl v 



A(p 

III ■ 

Pi) A(p 

in ■ 

x\) A(p 

mi %2 ) 

A(x 1 ,p 1 ) A(x 1 ,x 1 ) A(xi,x 2 ) dx 
A(x 2 ,pi) A(x 2 ,x 1 ) A(x 2 ,x 2 ) 

~ A(p 

in i Pi) A(p 

mi %l) 



+ 



4>(xi,x 2 ) det 



hxl v 



A(x 2 ,pi) A(x 2 ,x l ) 



dx. 



In view of Lemma 12.31 and (|2.2ip in Lemma 12.41 below we obtain (|2.19p . 



Lemma 2.4. Suppose that f : 



has a continuous derivative and that 



g : R — ► K has continuous first partial derivatives. Assume that 

|/(x)|, |/'(x)| < Ce- 

\g(x,y)\, \g'(x,y)\<Ce-*-y. 
Then, for 1 < i,j < m, it holds that 



(2.20) 

(2.21) 
Proof: 



1 (x-Pj) 2 

e 4, f(x)dx = f(pj) 



'Aire 



1 



iZEi \/47r 



e" 2 rdx + 0( v / e) 



>Ii J Ij V 47T6 



1 (g-y)j 

e 4e 



5(^,y)dxdy= / fi((x,x)dx + 0(v / elog< 



'47re 



e 4 E /(a-) dx 



Pj 



e " 4 /(Pj + Vez)dz 



By Taylors theorem 

/fo + V~ez) = f(Pj) + V~ezf( Pj + 9 e (z)) 
where 9 e (z) is a number between and \fe~z. Since by assumption 

\f'( Pj + e e (z))\ <ce-v+^w 

we obtain ([2T2TJ]) . 
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roo roo 



1 (x-y) 2 

e 4£ g(x, y) dxdy 



y-x 



Pi 



e 4 g(x, x + yfkz) dxdz 



By Taylors theorem 

g(x, x + \fez) = g(x, x) + \/ezg'(x, x + 9 e (x, z)) 
where 9 e (x, z) lies between and y/ez. The error can be discarded since 

.2 



/•oo /*oo 



e 4 |^</(x,x + 9 € (x,z))\ dxdz 



< C / dx / -—^e-V^+v^M dz < C. 

Jp t J-oo V47T 



We now split the main term into two terms. 



dx 



Pi 



dz 



xA 7 



e 4 g(x, x) 



4tt 



/•Pi-yeloge poo 

/ dx / 

poo poo 

/ dx L-* 



dz^=e * g(x,x) 

V47T 



+ 



We can estimate the first integral by 



dz -= e 4 g(x,x) =: / + 
V47T ./ 1 ./ 2 



./l ip 



dx 



eke™ -2 * < -C\/^loge. 



If x > pj — y^e log e then < C log e and hence 

roo j 2 



e 4 dz 



1 -£ , 

e 4 dz + 



47T 



J-oo \/4~7T 



(loge)' 

e T = 1 + O { e 4— 



We finally get 



Pi— v/i log e 



(l°g£) Z 

1 + e" 4 



x) dx 

/•oo 

= / g(x,x) dx + C(\/ilog< 



This concludes the proof of the lemma. 
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3. Theorem [O] 

3.1. Multi-layer discrete PNG. Before we give the proof of Theorem 1 1.31 
we must present some preliminary results. 

How does one get a hand on the process h described in the introduction? 
In [7j it is shown that h can be embedded as the top curve in a multi-layer 
process given by a family of non- intersecting paths {hi, < i < N}, h = ho. 
It turns out, see [7], that this multi-layer process is an example of a discrete 
determinantal process. 

Theorem 3.1 (Johansson). Let u,v S Z be such that \u\,\v\ < N and let 
q = a 2 . Set 

G(z w) = (l- a ^W^q-a^- 

and 

~ 1 f dz f dw z . 

K N (2u,x;2v,y) = / — / G{z,w) 

(2m) z J lT2 z J 1ti w z-w 

where j r is the circle with radius r centered around the origin, a < r\ < 
T2 < 1/ct and x,y 6 Z. Furthermore, define 

fouM^V) = ^tf j {y - x)6 G (e id ,e ie ) d9 

for u < v and 4>2u,2v(x, y) = for u > v. Set 

K N (2u,x;2v,y) = K N (2u, x; 2v, y) - (p 2 u,2v(x, y). 

Then, 

F[(2u,x] u ) £ {(2t,hi(2t,2N - 1));0 < i < N, \t\ < N} , 

\u\ < N, 1 < j < k u ] 
= det (K N (2u,x!--2v,xf )) MM<NA ^ A ^ 
for any x 2u G TL and any k u G {0, . . . , N}. 

The asymptotic information about the kernel needed to prove Theo- 
rem [L3] is contained in two lemmas. The first can be extracted from chapter 
four in [7] and the proof of the second is provided at the end of this section. 
Please note that we make a slight redefinition of the function (j) from the 
last section. However, for the purposes of this text <p acts as one and the 
same. 

Lemma 3.1. Let t,t' be any real numbers such that 

u = i±^ d -lAT2/3 r G Z+ 

I — a 

v = l±^ d -^ N vz T > € z 

1 — a 
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Let x, y G Z + and define x' , y' by 

x = 2a(l - a) -1 JV + (x' - T 2 )dN 1/3 
y = 2a(l - a^iV + (y' - T ,2 )dN 1/3 . 

For any L G R i/iere exisi positive constants, c and C , such that 

\K N {2u,x;2v,y)\ < CN~ 1 ^ 3 e~ c ^ x ' +y '^ 

ifx',y'>L. 

If \x'\, \y'\ < log TV, then the exists c > such that 

dN 1 ' 3 K N {2u, x; 2v, y) = e z V^ v -^ A(t, x; t', y') + 0{N' 

Lemma 3.2. Let x,y G Z + and define x',y' by 

x = 2a(l - a)" 1 ^ + x'diV 1/3 
y = 2a(l - a)- x iV + y'dN 1 / 3 . 

Take s > 0, let u ~ iV 2 / 3 and define v by 

f = u H a siv 7 

1 — a 

where < 7 < | . There exists a constant C > snc/i i/iai 



where 



and 



1 



(7/7/) = ; e 4siVT-2/3 



/or a// x, y. 



3.2. Proof of Theorem 11.31 This proof is really a discrete analog of the 
proof of Theorem 11.11 Unfortunately things are more involved in this case 
where TV 7-2 / 3 plays the role of e. 

Please recall that J\ = fiN + ipdN 1 / 3 where // = 2a (1 — a) -1 and q = a 2 . 
Set Ji = Jj_i + G Z, i = 2, . . . , m, and 

Ii = {z G Z|z > J;} . 

Here the y^'s are arbitrary numbers such that Ji G Z. For later convenience 
we also define ipi, i = 1, . . . , to, by Ji = fj,N + ipidN 1 / 3 . 
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We will prove that 

P [# -h = ■ ■ ■ = #Jrn = 1, #h = ■ ■ ■ = #ln = 0|#Ji = 1, = 

= <p2K u 2K 2 (Jli J?)'" (f>2K m -i,2K m (Jm-l, Jm) 

+ o([N-"il 2 ) m - 1 N- A j . 



This implies Theorem 11.31 

4>2K 1 ,2K 2 {-h, J2) ■ ■ ■ §2K m - X ,2K m {Jm-\, Jm) 

1 -A. 1 1 



-e 4s 2 e ^ ; (l + 0(N~ c )) 

by Lemma 13.21 The sum of this function over the sets Ai is a Riemann sum 
that is well approximated by the integral in Theorem 11.31 
Define the finite integer intervals Ij, 1 < i < m, by 

h = {z G Z; Ji < z < [/iN\ + N} . 
The probability of finding a particle in Ii but outside of Ii is very small: 

P \ ii) > l] < P[#x = l]= 

^J + (> 2/3 + ^) diVl/3 ) 

<Ce-^ 2/3 ^e-^73 =0 



fc=0 

This means that we can work with Ij instead of I{. We now proceed much 
like we did in the proof of Theorem 11.11 If we set 

A = {# Ji = 1, . . . , # J m = 1} , 

then 

PL4,#Il = • •• = #Im = 0] +P[A#/i = 0,(#J 2 = ... = #/ m = 0) c ] 

= P[A #/x = 0] 

where 

F[A, #h = 0, (#J 2 = . . . = #/ m = 0) c ] 

= P[A,#Ji = 0,U£ 2 {#I i ^0}] 



i=2 i=2 
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and 

P[A, #h + #h] = E[X{#J 1= i} • • • X{#J m =i} • X{#h^#h}} 
= E[# Ji • • • # J m ■ X{#h^#h}] 
<E[#J 1 --.#J m (#/ 1 -#/,) 2 ]. 

The second equality holds since the probability of finding two particles at 
the same place is zero. 

We need to prove three things: 

f. p[A#/i = o] 

= <hK 1 ,2K 2 (Jl, J2) ■ ■ ■ <hK m - U 2K m {Jm-l, Jm) 

xF[#J 1 = l,#I 1 = 0] 
+0 (N- 1 / 3 - c (N-~</ 2 ) m - 1 ^ 

2. E[# Ji ■ ■ • # J m (#7i - #Ii) 2 } = O (AT-V3-c (Ar -T/2 r -i^ 

3. P[# J x = 1 , #/! = 0] > CiV" 1/3 

Before giving the proofs we need some preliminaries. 

When summing a function f(x) over, say, I\ we can write 

where Ti ~ iV. The next lemma will be frequently used later on. 
Lemma 3.3. There exists constants C\,Ci > suc/t f/mi 

00 

^0^/^/3^^ N -l/3 < Ci 
fc=l 

and 

N 2 

fc=i 

for any x € M. 
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Proof. 



oo oo 



k=l 



'k-xN 1 ' 3 



.or 1 ' 3 



oo 

5 £ ■ 

fc = — OO 



ifc-zJV 1 / 3 



, N- 1 ^ = \f := xN 1 / 3 - [xN^ 3 \ 



< 



fc=— oo 

(j^ 75 ' ) w " 1/3+,A JV_1/3+ §' A JV "" 

fe=l 



fc=— oo 



JV 2 



c ATl/3_7V7 



J^jg (k/N l / 3 ,x\ < CN' 1 



1 



1 



fc=l k=l 

< CN^ 1 log iV + CN^ /2 + CJV~ 7 log JV < CN-~< /2 



□ 



We now turn to the proof of 1. As in the proof of Theorem ll.il we get 



x ' l) k 



(3.1) F[A, #h =0}=J2 -JT E * Jl • • • * J ™ # J i 



k=0 



rl*] 



For < r < m — 1 set 
D r (ife) 

= 4>2Ki,2K 2 (Jli ■h)4 l 2K2,2K3(J2, >h) ' ' ' <p2K r ,2K r +i(Jr , Jr+l) ^ 

Xi£li,l<i<k 



K(J r+ i,Ji) K(J r+ i,J, 



r+2) 



K{J m ,J\) K{J m ,J r+ 2) 
K{x i ,J 1 ) K(xi,J r+2 ) 



K(J r+ i,J m ) K(J r+ i,Xj 

K {Jrnj Jm) K ( J m , Xj ) 
Ki^XiiJrri) K{Xi^Xj) 



The indicies i,j run from 1 to A; and if r = the (empty) product of <j>- 
functions is to be interpreted as 1. Let D r (k) be like D r (k) but having 
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K(J r+ i, J r +2) in position (1,2) in the matrix. We want to show that 
\D (k) - D r (k)\ < AT-V3-C (V^ 2 )" 1 " 1 (Ck)^ 

which, by the induction argument in the proof of Theorem ll.il follows if we 
can prove that 

(3.2) \D r (k)\ < N- 1 ^ (N-^) m ~ l (Ck)^. 

To show this we shall use Hadamard's inequality and therefore need to es- 
timate sums of column elements squared (confer with the proof of Theorem 
1). Lemmas I3.1| 13.21 and 13.31 will be frequently used below. 
Column 1: 

m k 

K 2 {J u J 1 ) + Y,K 2 {x i ,J l )<CN~ 2 '' i {m + k) 

i=r+l i=l 

Column 2: 

m 

K 2 (J r+1 ,J r+2 ) + K 2 {J u J r+2 )<CN- 2 ^m 

i=r+2 

and 

k 

Y K2 ( x i> J r+2) 
i=l 

k 

< CN~ 2 ' 3 Y [l + <t> (k/dN 1 / 3 ,^ - t/v +2 ) 

i=l 

+ N l ^(t> E (li/dNV 3 ,lh - VV+ 2 ) ] 2 ■ 

Columns 3, . . . , m — r (r + 3 < j < m), if they exist: 

m k 

Y ^{JuJ^ + Y^i^Jj) < CN-~<{k + m) 

i=r+l i=l 

Last k columns (1 < j < k): 

m k 

K 2 {J i ,x 3 ) + Y J K 2 {x l ,x j ) < C(k + m)N- 2 / 3 e- cl > N ~ 1/3 

i=r+l i=l 
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Using Hadamard's inequality we get after some manipulations that 



\D r (k)\ 



Ti 



h,...,l k =l 

k 

X Z [l + 0(ii/^ 1/3 ,^-^r+a) 



8=1 



i=l 



It follows from Lemma 13.31 that 



Ti 



£ e- dl7V " 1/3 (li/dN 1 / 3 ,^ - ^ +2 ) JV-V3 < C 



and also that 

3i 



Ji=i 



<f> E (li/dN 1 / 3 , ^ - W+ 2 ) | < CiV^ 2 



From this we get (|3.2 j) . 

To get 1 we also need to show that 



(3.3) det 



K(J m ,J x ) K(J m ,Xj) 
K(x i ,J 1 ) K(xi,Xj) 



l<i,j<k 



E det 

Xidl\,l<i<k 



K(J x ,Ji) K(J liXj ) 
K(xi,J{) K(xi,Xj) 



l<i,j<k 



+ N-^- c O ((Cfc)^) 



Write 



diV 1 / 3 



and consider first the case 1 < h < A rl/3 log N. From Lemmas EE2 and O 
it is straight forward to deduce that if z = x% or z = J\ then 



K(J m , z) = K(J U z) + 0{N- l ^- c ). 
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We now expand the determinant in the sum to the left in (13.30 . 



det 



K(J m ,Jx) K(J m ,Xj 
K(xi,Ji) K(xi,Xj] 



det 



+ 0(7V-V3- C )^ de t [ K (xi, Ji) ^(^,x i )] 1 < iii < feJ 
P =i 

+ 0{N-V 3 - c ) det 
We now use Hadamard's inequality to get 

Ari/3 log jv 

|det [Kfa, Ji) K { x u x j)]i<ij<k,j^ P \ 



l<i,j<k 

K{J lt J x ) K{J x , X j) 
K(xi,Ji) K(xi,Xj) 



l<i,j<k 



'Jl\l<i,j<h 



< ^ (CkN~ 2/3 ) k/2 e~ N ~ 1/3( - h+ - +lp - 1+lp+1+ -- +lk * ) 

< (Ck) k/2 log N 



and 



Ari/3 logA r 

J2 |det[ J PC(x i ,x i ] 1 <. j .< fc 



ii=l 

This takes care of the summation over 1 < li < iV 1//3 log N, 1 < i < m. By 
using Hadamard's inequality once more one readily shows that the contribu- 
tion coming from the remaining terms in the sums in (|3.3[) is small enough 
to make (|3.3p hold. 

We now prove 2. Note that 

- #/ 4 ) 2 = if + 1? +h + h- 2Iih. 
By arguing as in the proof of 1 above we obtain 



E 
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<f>2Ki,2K 2 (Jl, h) ■ ■ ■ <l>2K m -i,2K m (Jm-l, Jm) 

K(J m ,Ji) K(J m ,x s 



x det 



Kr,s<fc 



K{x r ,J\) K{x r ,x s 
where u = l,i and k = 1,2. One also gets 

= 4>2K 1 ,2K 2 (Jl, Jz)-" 4>2K m - 1 ,2K m (Jm-l, Jm) 

K(J m ,Ji) K(J m ,x) K(J m ,y) 
x ^ det K(x,Ji) K(x,x) K(x,y) 
xeii,yeli [ K(y,J{) K(y,x) K(y,y) 

+ 0(^-1/3—^) 

We omit the details. Using Lemma 13. II and Lemma 12.31 one readily gets 



E det 



K(J m ,Ji) K(J m ,x s ) 
K(x r ,Ji) K(x r ,x s ) 



Kr,s<k 



E 



#Ji#/f ] + 0(7V 



-1/3-CN 



for k = 1,2 and 

E det 



K(J m ,Ji) K(J m ,x) K(J m ,y) 
K(x,Jt) K{x,x) K(x,y) 
K{y,J{) K{y,x) K{y,y) 



E 



#j 1 #4 2] +e>(iv- 



-1/3-0 



We now see that 2 follows if 
(3.4) ^2 4>2K 1 ,2K l (x,y)det 



K{J m ,Ji) K(J m ,x) 
K{y,J 1 ) K(y,x) 



= E[#J 1 #/ 1 ] + 0(iV-V3- C) . 

We shall prove this by showing that both sides are well approximated by 
integrals. On the integral containing the function <p we can then apply 
Lemma 12.41 

By using Lemma 13.31 we get rid of the error term associated with <j)E'- 



N 



E M^ 1 + ^T73'^ + 



h,h=l 



e n 



73 2/3 



A/ 



< C £ e"^ AT-2/3-7/2 < Cjv -l/3- 7 /2 
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The following calculation, again using Lemma 13.31 shows that the main 
contribution to the sums in (13. 4j) comes from summing over 1 < h,fo < 
ArV3 log AT. 

N N h+h 

; 1= jVi/3i og jvZ 2 =l 

< jh Ce- h / Nl/3 N-^ < ON' 1 

We shall use Euler's summation formula for two variables: 

Lemma 3.4. Let f(x, y) be a function of two variables such that its partial 
derivatives up to second order are continuous in the rectangle 

{(x, y)\a < x < b, c < y < d} 

where a, b, c, d are integers. Then 

rb r d 



a<m<b c<n<d 



b rd 



po pa 

Y Y /( m ' n ) = / / /(>>y) dxd y 

™ J, ^ J " Q J C 

I 

fx{x,y)(x - [x\)dxdy 

f 

+ / / fy{x,y)(y - LyJ)dxdy 



a J c 
b r d 



a J c 
b ,-d 



+ / / fxy(x,y)(x - [x\)(y - |j/J)dxdy 



i a J c 

The case that we are interested in is when 



f(x, y) = <f> (Vx + V* + ^) g (x/dNWy/dN 1 ") iV" 1 

where 

\g x (x,y)\,\g y {x,y)\,\g x {x,y)\ < Ce~< x+ ^ . 

We need to show that the integrals involving the absolute values of f x (x, y), 
f y (x, y) and f xy (x, y) are negligible. We only present the details for \f x (x, y)\ 
here, the other terms are treated similarly. 

N 1 /^ log N pN 1 /^ log N 

\fx(x,y)\ dxdy 

< (dN 1 / 3 ) 2 / / f x ({x - ^dN 1 / 3 , (y - iP^dN 1 / 3 ) dxdy 

POO f'OO 

< CN- 2 ' 3 / / (\Mx,y)\ + <Kx,y)) e~ c ^ dxdy 
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By Lemma ([27 



oo roo 



<j>{x t y)e~ c{x+y) dxdy < C. 



lip! Jlpi 

The remaining term demands some analysis. 



\Mx,y)\e- c(x+y) dxdy 



oo r-co 



\x - y\ 



oo 

dx 

i/'i 



i 27V7-2/3 V4vriV7-2/3 

1 



^ (x-y) 2 

g 4JV 7-2/3 g-c(a!+») dxdy 



+ 



2iVT-2/3 V4^iV7-2/3 



e 4Jv7 -2/3 e -c(x+y) dy 



dx 



x-y 



^ 2Af7-2/3 V4vriV7-2/3 



e 4JVT-2/3 e -c(x+J/) dy 



dx 



01 



</>{x, y)e~< x+y ^ +c </>(x, y)e-< x+y ^ dy 

< CN 1/3 ^ /2 + C < CN 1/3 ~ l/2 



We can do the same calculation for the remaining integral. The \ f x (x, y)\ 
integral is hence O (jV -1 / 3-7 / 2 ) and the same goes for the \f y (x,y)\ and 
\fxy(x,y)\ integrals. 
Set 

A T1 (x, y) = A(r l ,x + rf; n, y + r 2 ) . 

Applying the above calculations to the left hand side of (|3.4p and using 
Lemmas E21 E2Q and O we obtain 



K(J m ,Ji) K(J m ,x) 
K{y,J t ) K(y,x) 



^2 4>2K 1 ,2K i (x,y)det 

xeh,yeh 

N l/3lo S N 



1 



«1,«2 = 1 

e~ T1 ^^ A^^ + h/dN 1 ' 3 ,^) A T1 (^i + h/dN 1 / 3 ,V>i + h/dN 1 / 3 ) 

+ 0(N- 1 / 3 - c ) 



1 



d/VV3 



det 



i/'i 



A T1 (x,^i) A T1 (x,x) 



dx + C(iV 



-1/3-cx 



We get the same expression for the right hand side of (|3.4j) when applying 
Euler's summation formula. This concludes the proof of 2. 

Let F2(t) be the Tracy- Widom distribution function corresponding to the 
largest eigenvalue of the Gaussian Unitary Ensemble (GUE), [15]. That 3 is 
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true follows from the fact that F^t) > Vi, see [IS], together with the next 
lemma. 

Lemma 3.5. Let J\ and I± be as above. It holds that 

P[# Ji = 0, #A = 0] = ^73^2(^1 + rl) + 0(iV-2/3). 

Proof: This will, again, be an exercise in using Hadamard's inequality. We 
have the following representation for F' 2 (see the third equality in (|2.14p ): 



(3.5) i? 2W=y"^T^/ det(^(x i ,x j ))o<i, j < fc d fc x 

where xq = t. In three steps we will now show that 
00 (—i) k 

dNl/3 2 det(K(x l ,x j )) <i, j <k 

k=0 ' a; 4 eJi,l<i<fc 

where xo = J\, is well approximated by the right hand side in (|3.5|) . By 
(|3.ip this will prove the lemma. In steps one and two we will use Lemma 
13. II to insert the kernel A instead of K. In in the last step we show that we 
can change from summation to integration. 

First we show that we can sum over xs = fj,N + (^1 + U/dN 1 / 3 ) dN 1 / 3 
where 1 < Zj < iV^logiV, 1 < i < k, instead of over I\. By Hadamard's 
inequality and Lemma IBTTl 

det(^(x i ,x i )) <. J .< fc < YlJ2 K2 ( x i' x i< 

\j=0 i=0 

1 1 2 

< ( C(k + l)N- 2 ' 3 \{C{k + i) N -V* e -ch/N^ 1 

k 



< N- i ' 3 (c{k+ 1))^ n e_dj/ ^ 

j=0 



We have that 



00 k TV 1 / 3 log TV k 

Y Yle^/^N- 1 / 3 - Yj Yl e ~ clj/Nl/3N ~ 1/3 

l<i<fc l<i<k 

ATi/3 logA r ^ k 

^ k E J2H e ~ ch/Nl/ * N ~ 1/3 ^ kCkN ~ 1 - 

2<i<k 
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Since 

fc=l 

we see that we can indeed restrict the summation. 

In the second step we replace K by A. As before we shall use the notation 
A T (x,y) = A(x + r 2 ,y + r 2 ). For 1 < k < N 1/3 logN it holds by Lemma 
GO that 

det(K(xi,Xj)) <ij<k 

- m^ det ( AT ' (k/dNl ' 3 ' ll/dNl,3) + ° (N ^ 

where we let Iq = ipidN 1 / 3 . If we expand the determinant in the right hand 
side we get (k + l) 2 error terms of type 

AT — C , . 

— ^^det (A^k/dNWydNW) + 0(N-*)) ^ k . 

An application of Hadamard's inequality together with Lemma 13.11 shows 
that the total error we get when changing from K to A T1 is of order ]V — V 3- c . 
We omit the details. 

Finally we want to go from summation to integration. To do this we shall 
use that 

JVVsiogJV 

(3.6) Yl A^(l i /dN 1 / 3 ,x)A^(y,l i /dN 1 / 3 ) 



/•oo 

dN 1 ' 3 / A T1 (z,x)A T1 ( y y,z)dz + (e~ x - y ) 
Jo 



and 

JV x / 3 logiV 

(3.7) Yl A^{k/dN l / 3 ,k/dN 1 / 3 ) 



/•oo 

dN 1 ' 3 / A Tl (z,z)dz + 0(l). 
Jo 



This follows from Euler-Maclaurins summation formula and Lemma [2. 31 We 
will show that 

ATi/3logAf 



l, 

Ki<k 



' (di\Tl/3)fe+l V VV ' J ' J J0<i,j<k 

+ ({Ck) 1 ? N- 2 ' 3 
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where Iq = dN 1 / 3 ipi and yo = V'l- This will prove the lemma since 
For r = 0, . . . , k we set 



k=l 



oo. 



D, 



(dN 



./3)*-r+l det (^ Tl ^'))o<iJ<fc 



where 

r Vi i = 

= < j/i 1 < i < r 

{ k/dN 1 / 3 r+l<i<k 

Please note that Do is what we sum over in (I3.8P and that is what we 
integrate over. D r should roughly be what we get after having changed 
summation over 1%, . . . , l r to integration over y%, . . . , y r . We can expand D r 
in such a way that we get k 2 terms of type 

± (djvl/ 3 )fc _ r+1 ^ ri (^cWl/^V 71 (Ir+l/dN 1 / 3 , Z jo ) 

x det (A Tl (zi,Zj))o<i t j< k 

j¥=r+l,3o 

and one term 

rl/3 )fc -r+l AT1 (Ir+i/dN 1 / 3 , l r+1 /dN^ 3 )det (A Tl (z h Zj ))o^ 



(dN 1 



<k ■ 



We now apply f|3.6|) and (|3.7p and therefore need to deal with the corre- 
sponding errors. 



( 



3¥=3o,r+l 



1/2 



< C(N 1/3 ) k ~ r+1 e~ Zl o~ z J- 



. 3=0 

VVio,r+l 



/ 



<(7V 1 / 3 )*- r+1 (C7(fc-l))^ JJ e 



3=1 



Since 



JV 1 / 3 log N k 



j d r x j2 n e ~ czj ^ c k (N i / s ) k -^ 



l i= l j=l 

r+2<i<k 



we find that the error from the k 2 terms of the first type is estimated by 

k 2 (C{k-l)) k -^N- 2 ' 3 . 
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The error coming from the remaining term can be treated in the same way. 
Changing from summation over li to integration over j/j, 1 < i < k, hence 
results in an error estimated by 

kk 2 (C(k - l))^iV-2/3 = ( C Q±p N -2/3 

as needed. 

Proof of Lemma 13. 2t By definition 



(1 - a) 2 ("- u ) f 77 
2tt 



ji(y—x)9+(u—v) log (1+q 2 — 2a cos 9) 



4>2u,2v{x,y) 

Define 

g{6) = log(l + a 2 - 2a cos 9) 

in [— 7T, 7r]. This function is analytic in a neighbourhood of zero and a 
Maclaurin expansion gives 



g(6) = log(l-a) 2 + 



a 



T e 2 + c 2 # 4 + o(# 6 ) 



(l-a) 2 

where C4 < 0. It is easy to see that for any 5 > there exists e > such 
that 

g{6) > log (1 - a) 2 + e 

if 101 > 5. Hence 



(l-a) 



2(u-u) 



|0|><5 



2vr 



< 



1 

2^ 



We expect that the main contribution to </>2« 2« will be 



2?r 



2tt 



1 



1 

27rVi(iiV 1/3 
1 



e 



00 iV^sLt^-U 
e 



1 



(x'-y 1 ) 2 



diVV3 V4vrsiVT- 2 / 3 
Below we will analyze the error. For simplicity we take s = 1. 
Define fc(0) by 



This means that 



g{9) = log (1 - a) 2 + ^—^ (9 2 + . 

00 

h(e) = J2h k e k 

k=A 
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where /14 < 0. Note that h is even since g is and also that, for 5 small 
enough, h(6) < if \0\ < 5. The error becomes 



Err 



e i{y'-*')dNVH F{e)de 



where 



_ e -d 2 Nie 2 _ e -d 2 N-<6 2 ~d 2 N-ih(6) 



Next we integrate by parts. 



Err < 



i(y> - x^dN 1 / 3 



+ 



i(j^-x')rfjv 1/3 e 



F{6) 



< 



\y> -x'ldN 1 / 3 

3 -d 2 N"'S 2 



j & Sy'-z')dN^e F , {e)de 



\y' - x'ldN 1 / 3 



+ 



y> - x'\dN 1 / 3 7_a 



\F'{6)\de 



The last integral will be easy to compute if we can find out where F'{6) 
changes sign. 

F'(9) = 2d 2 N^9e- d2N ''^ +h ^ (l + ^ - e^^^ 
A point in [—6, 5] \ {0} where F' changes sign will satisfy 



h{9) 



d 2 N^ 



log 



1 + 



h'(0) 
28 



This shows that if N is large then F' has two zeros ±60 in [—5,5] \ {0}. 
Moreover, 9q is of order N~"'/ 2 . Given this information we check which sign 
F' has in different intervals and get 



\F'{0)\dO = 2 f \F'{6)\d6 
Jo 

= - / F'{6)d6+ / F'{6)d9 
Jo Je a 

= F(0) - F(0 ) + F(S) - F(0 ) 
= O(AT-T). 

This almost finishes the proof of the second inequality in the lemma. We 
should not forget the exponentially small error terms that appeared above. 
They do not have the factor \x' — y'l^ 1 in front of them. However, a couple 
partial integrations can be used to take care of this obstacle. 
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The first inequality in the lemma follows from the following calculation. 

/ \F(0)\dO= [9 = tN~^ 2 ] 
Jo 

dt 



ply 11 

= N^ /2 J e -d 2 t 2 -d 2 h(tN->/ 2 ) ^ _ e cPN-'h{tN-^ 2 )^ 

< N~~</ 2 / e-^ 2 (l - e- C2N " H ') dt 

< N~^ 2 J te-^ 2 (l - e - C2N "' t4 ) dt + 0(N 



We now use partial integration. 

N"t/ 2 8 



L 



te-^ 2 (l - e - C2N ~ lli ) dt 

N~t' 2 & 



e c,t 2 l _ e -c 2 N~^ 



1 

2ci 

C\ Jl 

= o (n-^' 2 ) 

This concludes the calculations in this section as well as in this paper. 
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